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Abstract 

We address, within the dual superconductivity model for color confinement, 
the question whether the Yang-Mills vacuum behaves as a superconductor of 
type I or type II. In order to do that we compare, for the theory with gauge 



{Nj ■ group 577(2), the determination of the field penetration depth A with that of 

^* . the superconductor correlation length £. The latter is obtained by measuring 

the temporal correlator of a disorder parameter developed by the Pisa group 
to detect dual superconductivity. The comparison places the vacuum close 
^ ■ to the border between type I and type II and marginally on the type II side. 

\& • We also check our results against the study of directly measurable effects 

such as the interaction between two parallel flux tubes, obtaining consistent 
indications for a weak repulsive behaviour. Future strategies to improve our 
investigation are discussed. 

a.; 

I. INTRODUCTION 

Color confinement emerges as a fundamental property of strongly interacting matter 
from experimental facts, like for instance the absence of fractionally charged particles. Even 
if lattice simulations provide evidence that confinement is realized in the theory of strong 
interactions, a full theoretical explanation of it starting from QCD first principles is still 
lacking. However models exist which relate confinement to some property of the funda- 
mental state of the theory. One of those models is based on dual superconductivity of the 
QCD vacuum [1-3]: according to this model confinement of color is due to the sponta- 
neous breaking of a magnetic simmetry which yields a nonvanishing magnetically charged 
Higgs condensate. The dual Meissner effect compels the electric field between static colored 
charges in narrow flux tubes, giving rise to a linearly rising potential and to confinement. 
The broken magnetic group is chosen by a procedure known as Abelian projection [4]: a local 
operator 4>(x) transforming in the adjoint representation is diagonalized, leaving a residual 
U(l) Nc ~ l gauge symmetry. 
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A superconductor is characterized by two fundamental parameters, the correlation length 
£ of the Higgs condensate and the field penetration depth A: they determine whether the 
superconductor is of type I (£ > A) or type II (£ < A). In a superconductor of type I an 
external field B is always expelled from the medium till a critical value B c beyond which 
superconductivity disappears. In a superconductor of type II there are instead two different 
critical values B c \ and B C 2, and for B c \ < B < B C 2 the external field can penetrate the 
medium in the form of Abrikosov flux tubes, without disrupting superconductivity. Another 
relevant property of type II superconductors is the repulsive interaction between two parallel 
flux tubes, which is instead attractive for type I superconductors. 

In the framework of the dual superconductor model, understanding whether the QCD 
vacuum behaves as a type I or a type II superconductor is an issue which can help clarifying 
the dynamics of color confinement and of flux tube interactions. The question can in principle 
be answered by QCD numerical lattice simulations and several efforts have been done in 
the past in that direction, mostly for the pure gauge theory with 2 colors. A direct way 
to determine A is a lattice analysis [5-10] of the flux tube which is formed between two 
static color charges: the longitudinal (chromo)electric field E z decays asymptotically as 
E z = AKo(j) at a radial distance d from the tube axis. 

The determination of £ is less straightforward: this parameter has been found in literature 
mostly either through an analysis of violations of the E z = AKq(j) behavior close to the 
center of the flux tube [6,10] or through some global fit to the whole set of Ginzburg-Landau 
equations [7-9]; a determination based on a direct analysis of the condensate distribution 
around the flux tube has also appeared recently [11]. An approximate picture has emerged 
placing the 577(2) Yang-Mills vacuum roughly at the boundary between a type I and a type 
II dual superconductor. 

In the present study we consider the case of 577(2) pure gauge theory and follow a 
different strategy, aimed at determining the mass of the Higgs field m# = l/£ through the 
analysis of the temporal correlator of an observable directly coupled to it: that is the operator 
\i developed by the Pisa group which creates a magnetic monopole (see Ref. [12] for a detailed 
discussion about its definition and also Ref. [13] and [14] for related parameters). Its vacuum 
expectation value (v.e.v.) (/x) is a good disorder parameter detecting dual superconductivity 
((/i) 7^ 0) and the transition to the deconfined - normal conducting phase ((//) = 0) both 
in pure gauge theory [15-17] and in full QCD [18,19]. We will compare results obtained for 
£ in this way with those obtained for A through the usual analysis of the field inside the 
flux tube. As a further independent method to characterize the QCD vacuum, we will also 
directly study the interaction between flux tubes by measuring the electric field in presence 
of two couples of static charges. Preliminary results concerning the determination of £ have 
been reported in Ref. [20]. 

In Section II we will review the definition of the disorder parameter (fi) and present 
a determination of £ based on the measurement of its temporal correlator. The results 
obtained for £ will be compared in Section III with those obtained for A. Conclusions 
concerning the typology of the vacuum will then be checked against a direct analysis of flux 
tube interactions in Section IV. Finally in Section V we will present our conclusions and 
discuss possible improvements as well as possible future extensions of our study. 
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II. DISORDER PARAMETER FOR DUAL SUPERCONDUCTIVITY AND 
DETERMINATION OF THE CORRELATION LENGTH f 



A. The disorder parameter and its temporal correlator 

A disorder parameter detecting dual superconductivity can be constructed in terms of 
an operator \x which creates a magnetic charge. It can be defined in the continuum as [12]: 

fi a (S,t) = exp( i JdyTi{<P\y,t)E{y,t)}b L {y-x)j (2.1) 

where 4> a (y, t) is the adjoint field defining the abelian projection, is the field of a monopole 
sitting at x and E(y, t) is the chromoelectric field. The construction of \x is analogous to 
that of a translational operator in quantum mechanics: it creates a magnetic monopole by 
shifting the quantum vector potential field by the classical field b±. On the lattice correlation 
functions of n(x,t) can be written as (see [12,15-17] for details): 

mmt^) = § = \ { ^ e ~Z ( 2 - 2 ) 

where S is the usual pure gauge action and S differs from S only at time slices t and t' . In 
particular, in the abelian projected gauge the temporal plaquettes 

n i0 (y,2/o) = US,yo)U (y + i,yo)U}(y,y + 6)U^y,y ) (2.3) 

are changed by substituting 

Wy, y ) - Ui(y, y ) = U^y, y )e iT6 i^ (2.4) 

where T is the diagonal gauge group generator corresponding to the monopole species chosen 
(T 3 = A 3 /2 is the only possible choice for the SU(2) gauge group) and b\ is the transverse 
vector field corresponding to the monopole (antimonopole) sitting at t (t r ) and x. 

The numerical study of the temporal correlator of \x as a mean to determine the monopole 
mass has already been considered for the U(l) pure gauge theory in 4 dimensions [12,21]. 
In the confined phase, where dual superconductivity is at work, ^ 0. Therefore at 
large temporal distances the correlator (fi(t,x)ii(0,x)) is dominated, by cluster property, 
by a term (fi) 2 plus a function which vanishes exponentially according to the mass M of 
the lightest state coupled to \i. Taking into account that we are computing a point-point 
correlator instead of a zero momentum one and neglecting the possible presence of excited 
states, we will consider as the simplest possible ansatz the leading large distance behaviour 
of the two point correlation function: 

e -Mt 

(fi(t, x )/i(0, x)) ~ (^) 2 + 7-^ • (2.5) 

Since the ratio of partition functions in Eq. (2.2) is an exponentially noisy quantity, it is 
not easy to measure the correlator directly and one usually measures: 
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p = ±\n(w) = (S) s -(S)s (2.6) 

where the subscript indicates the action that is used in the Boltzmann weight. The be- 
haviour expected for p at large t can be easily derived from Eq. (2.5); after introducing the 
adimensional lattice quantities M = aM, t = t/a, n = x/a and after rescaling 7 — > 027, 
where a is the lattice spacing, one obtains: 

pit) = — In /i t, ft )fi 0, n ~ — ^ 7 2.7 

ap (/x) 2 + je~ Mt / t 3 ' 2 

where 

Eq. (2.7) will be the basis for our fits to the temporal correlator pit), which will be dis- 
cussed in Section II B. Results obtained through a different observable, also related to the 
temporal correlator in Eq. (2.5) and introduced in Ref. [21], will be presented and discussed 
in Section II D. 

As a result of our fits we will obtain an estimate of £ M = aM" 1 . The fact that (p) is 
a good disorder parameter for dual superconductivity means that it is surely coupled to 
the condensing Higgs field. The natural expectation is therefore that £ M = £, which is true 
apart from the unlikely case where the actual field which condenses in the vacuum does not 
coincide with the lowest mass state having the same quantum numbers (in that case one 
would have £ < 



B. Monte Carlo simulations and discussion of results 

We have measured the correlator pit) using a magnetic charge defined in the so-called 
random abelian projection, which was proposed in Ref. [17] and is a sort of average over 
all possible abelian projections: in that case one thus does not need to perform any gauge 
fixing at all, with a great benefit in computational cost. The dependence of our results 
on the abelian projection chosen will be discussed in Section II C, where we will make a 
comparison with results obtained by taking the abelian projection in the gauge where the 
Polyakov loop is diagonal. 

The correlator pit) is composed of two terms (see Eq. (2.6)): 

pit) = (S) s - (S(t))s({) ; (2-9) 

since the first term is independent of t, we have only determined the expectation value of 
the modified action (S(t))§,^ : we notice that a different Monte Carlo simulation is required 
for each value of t. 

We have performed simulations at four different values of the inverse bare coupling, 
(3 = 2.4, 2.5115, 2.6, 2.7, in order to eventually check the correct scaling of our results to 
the continuum limit. For the determination of the physical scale we make reference to 
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the non-perturbative computation of the (3— function and to the determination of T c j yfd 
reported in Refs. [22,23], from which we have inferred the following values of the lattice 
spacing: a((3 = 2.4) ~ 0.118 fm, a((3 = 2.5115) ~ 0.083 fm, a((3 = 2.6) ~ 0.062 fm and 
a{(3 = 2.7) ~ 0.046 fm. The lattice volumes A 3 x N t have been chosen so as to have 
approximately equal spatial volumes at the three lowest coupling values: 12 3 x 16 at (3 = 2.4, 
16 3 x 20 at = 2.5115 20 3 x 20 at (3 = 2.6. At = 2.7 we have been compelled by 
computational constraints to use again a 20 3 x 20 lattice, which is smaller in physical units, 
but is however comparable, as for the spatial size, to a 12 3 x 20 lattice at (3 — 2.5115, where 
we have checked that finite size effects in the determination of £ are negligible, at least 
within our statistical uncertainties. Different values of the magnetic charge Q carried by the 
monopole have been used in some cases, in order to check that our results are independent 
of this quantity. 

The signal obtained for (S (£))§, & is mostly made up of a constant background: it is 
therefore essential to reduce the noise as much as possible to obtain a good definition of the 
exponentially decaying signal. In order to do that we have integrated analytically over the 
probability measure of each gauge link (as the other links were left fixed), thus obtaining 
an improved estimate for the local action density. The typical number of measurements 
taken for each determination of the temporal correlator has ranged from 10 6 to about 5 • 10 6 . 
We report in Fig. 1 a summary of the results obtained for the modified action density 

n(t) = (S(t))s {i) /6V. 

The expected behaviour for II(t) stems from Eq. (2.7) by simply adding a constant term; 
as a matter of fact, due to the high number of parameters in Eq. (2.7) and to the poor quality 
of our signal, we have been able to fit only the leading large t behaviour of Eq. (2.7),^ which 
taking into account the periodic boundary conditions in the time direction is 

(2.10) 

In Table I we report the fit results obtained for £ according to Eq. (2.10) as a function of 
the initial fitting point t ; since we are taking into account only the leading large t behaviour, 
and also in order to avoid contaminations from higher excited states coupled to /i, we must 
search for a plateau in £ as a function of t . That is not an easy task since, as we will soon 
discuss, the correlation length £ comes out to be of the order of 0.1 fm: that, combined 
with the very low signal/noise ratio characterizing our observable, leads to a signal which 
disappears after a few lattice spacings, so that fits for to > 4 are hardly feasible, except for 
the highest values of f3 and Q where £ is larger and the signal sharper. 

As a general rule, we have considered for our determination of £ the value of to after which 
the signal does not change considerably within errors: that corresponds to to = 2 for (3 = 2.4, 
to = 3 for (3 = 2.5115 and (3 = 2.6, t = 4 for (3 = 2.7. The corresponding values obtained for 
the x 2 test are the following: x 2 /d.o.f.(12 3 x 16, /3 = 2.4) = 0.46/3, x 2 /d.o.f.(12 3 x 20,(3 = 
2.5115) = 1.2/2, x 2 /d.o.f.(16 3 x 20,/3 = 2.5115) = 2.1/5, x 2 /d.o.f.(20 3 x 20,(3 = 2.6, Q = 



tA fit which takes into account also the next to leading term, e */£/t 3 / 2 , gives compatible results, 
but errors are of the same order of the fitted values. 



nffi = ^ + ^ W + - — rr/2 
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2) = 2.7/3, x 2 /d-o.f.(20 3 x 20,(3 = 2.6, Q = 8) = 1.3/3 and x 2 /d.o.f.(20 3 x 20,(3 = 2.7) = 
6.4/4. The two determinations obtained at (3 = 2.5115 on the two different lattice sizes are 
in agreement, thus indicating that finite size effects are not important within our present 
statistical errors. The two determinations obtained at (3 = 2.6 for two different values of the 
monopole charge are nicely compatible within errors, thus showing no significant dependence 
of £ on Q. 

For the two lowest values of (3 the resulting value of £ is quite close, actually compatible, 
with the value of the lattice spacing itself, so that lattice artifacts could play an important 
role. The situation improves for (3 = 2.6 and (3 = 2.7, which can then be considered as more 
reliable determinations: we notice that the two values to = 3 and to = 4 used respectively 
for the determinations at (3 = 2.6 and (3 = 2.7 are approximately equal when converted in 
physical units. 

In Table II we report a summary of the values in physical units obtained for £ as a 
function of (3, together with the lattice spacing a{(3). Our results seem compatible, within 
errors, with the correct scaling to the continuum limit; however, taking into account that 
the determinations at the two lowest values of (3 are very close to the ultraviolet cutoff 
and that in general our statistical uncertainties are still large, a reliable extrapolation to 
the continuum limit is still not possible. Rather we give as our best determination of the 
correlation length £ that obtained at (3 = 2.7, which, using a conservative estimate for the 
error, is £ = 0.11 ± 0.02 fm. We will come back to these results in Section III where we will 
compare them with those obtained for the dual penetration length A. 

C. Independence of the abelian projection 

In the present Section we will discuss the possible dependence of £ on the abelian projec- 
tion chosen to define /i. The natural physical expectation is that £ be an universal quantity 
characterizing the Yang-Mills vacuum, hence independent of the particular abelian projec- 
tion chosen. This is consistent with 't Hooft ansatz that all abelian projections are equivalent 
to each other: that equivalence also emerges from numerical determinations of (//), which 
have clearly showed that {fi) being zero or non zero is a gauge independent statement [15-17]. 
A possible theoretical argument is the following: the operator [i defined in one particular 
abelian projection creates a magnetic charge in every other abelian projection [19,24]; this 
implies that the lowest mass state coupled to fi should be universal, i.e. £ should be inde- 
pendent of the abelian projection chosen. In order to test that hypothesis we have repeated 
our measurements for the abelian projection defined by diagonalizing P(n, t) on each lattice 
site, where P{n, t) is the Polyakov loop at the spatial site n starting at time t. 

The updating procedure in this case is not as simple as in the case of the random abelian 
projection: changes in Polyakov loops modify the abelian projection and as a consequence 
also the modified action S, which therefore is not a linear function of the temporal links. 
On those links usual heat-bath or over-relaxation updatings are not possible and we have 
used a metropolis algorithm. Numerical strategies for noise reduction like link integration 
are no more feasible and as a consequence there is a considerable increase in computational 
effort with respect to the case of the random abelian projection. 

We have performed a numerical simulation in the Polyakov gauge at (3 = 2.4 on a 12 3 x 20 
lattice. The results obtained in this case show a good agreement with the determination in 
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the random gauge, as can be appreciated from Fig. 2. A fit according to Eq. (2.10) gives 
= 2.4) = 1.3 ± 0.8, which is in agreement, even if within the large errors, with the value 
obtained in the random gauge. 



D. Comparison with a different approach 

An alternative way to study the temporal correlator of the disorder parameter, inspired 
by studies in gaugeball spectroscopy [25], has been introduced for the U(l) pure gauge 
theory in Ref. [21], and consists in considering a new observable p, which is the derivative 
of \n{p(t, n)p(0, n)) with respect to the adimensional temporal distance t in place of the 
inverse gauge coupling [3. The expected behaviour for p can be easily derived from Eq. (2.5) 

d / - 3 \ ie~ m /i 3/2 

m s j-HMUMo,*)) - - (m + -.) W2 ; 7 ; iVt/ . 3/2 . (2.1D 

Two features of the new observable p are apparent from Eq. (2.11): 1) the disconnected 
large distance contribution has disappeared in taking the derivative, so that p takes into 
account only the interesting connected piece without any noisy background; 2) there are 
only 3 parameters (M, 7 and (p) 2 ) to be fitted. Both things could contribute to make p a 
better observable than p in order to extract the correlation length £ = M~ x . To test that 
possibility and to have an independent check of the results presented in Section II B, we 
have repeated our determination of £ by measuring the correlator p{t). 

A drawback of p is that its definition on the lattice requires a prescription for the dis- 
cretized derivative, implying the possible presence of further systematic effects due to the 
finite lattice spacing. Let us rewrite the (pp) correlator in the form 

<M%(0)> = H s \ vu)e . ps , (2-12) 

where ASq and ASf indicate the changes in the action in correspondence of the monopole 
creation and destruction operators respectively; we will define p by taking the symmetric 
derivative, which can be written as (see Ref. [21] for details): 

tii) = -l(<,M l+1 -AS t _ 1 )) s+ASo+ASi (2.13) 

We have performed our measurement at /3 — 2.5115 on a lattice 12 3 x 16. Our results are 
reported in Fig. 3 together with the best fit result obtained using Eq. (2.11) after taking 
into account the periodic boundary conditions in the time direction*: the first point included 
in the fit has been t = 3. We have obtained £ = 1.32(25) (£ = 0.110(21) fm in physical 
units) with x 2 /d.o.f. = 1.8/3. We conclude that the agreement with the results obtained 
by measuring p is very good (see Tables I and II): this consistency gives an indication that 



* Since p(i) is the first derivative of the temporal correlator, it is an odd function with respect to 
t = Nt/2, as is clearly verified from Fig. 3. 
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systematic effects in the determination of £ are under control both when p or p are used as 
observables. The precision on £ is similar to that obtained in Section II B: since a comparable 
statistics has been used, we conclude that the benefit of dealing with a connected observable 
is not very significant. 

III. DETERMINATION OF A AND TYPOLOGY OF THE VACUUM 

Several consistent determinations of the parameter A can be found in the literature [5-10]: 
in this Section we will present an independent one obtained from the study of the flux tube 
profile between two static color charges. 

We have analyzed, for two different values of the inverse gauge coupling, (3 = 2.5115 and 
(3 = 2.6, the abelian projected flux tube formed between a quark and an antiquark placed at 
16 lattice spacings apart from each other (corresponding respectively to 1.33 and 0.99 fm); 
the Maximal Abelian gauge has been chosen to define the abelian projection. In particular 
we have studied the correlation of the plaquette operator with a Wilson loop W(R,T) with 
R = 16 and T = 6, on a lattice 24 x 24 x 32 x 24, with the longer dimension along the qq 
axis. We adopt the following prescription for the electric field [10]: 



which is the definition satisfying the Maxwell equations on the lattice [10]. The flux tube 
profile has been studied at half a way between the two charges and at an equal temporal 
distance between the creation and annihilation times of the qq pair. The Jackknife method 
for correlated quantities [26] has been used in the statistical analysis. 

Having chosen a quite long flux tube, noise reduction is a critical point of our computa- 
tion: we have adopted a standard cooling procedure, looking for a stable plateau in A as a 
function of the cooling steps performed. 

We have fitted our data according to the solution E z = AKo(i) of London equation 
V 2 E Z = tsE z : since that is expected to be valid beyond a certain distance d from the flux 
tube axis, where the effects of the non superconductive core are absent, we look for a plateau 
of A with respect to the minimum distance c?o included in the fit. 

In figure 5 we report the dependence of A measured at (3 — 2.6 both as a function of the 
fit starting point do at a fixed number of cooling steps (N coo i = 6) and as a function of the 
number of cooling steps at a fixed fit starting point (do = 3). A plateau is visible in both 
cases and we choose d = 3 and 6 cooling steps as a reference. 

Our fitted value at (3 = 2.6 (Fig. 4, left) is X((3 = 2.6) = 2.58 ± 0.12; a similar analysis 
at (3 = 2.5115 (Fig. 4, right) leads to \((3 = 2.5115) = 1.96 ± 0.08. Converting our results 
into physical units (see Section II B) we obtain A = 0.163 ± 0.007 fm at (3 = 2.5115 and 
A = 0.160 ± 0.007 fm at (3 = 2.6 lattice in good agreement with previous literature (as one 
of the latest determinations we report A = 0.157 ± 0.003 from Ref. [10]). 

In Fig. 6 we report a summary of the results obtained for A and £ at the different 
values of the lattice spacing. While apparently £ is consistently lower than A, it is anyway 
clear that the two quantities are comparable, in agreement with the findings of previous 




(3.1) 



(ti(W AbPr (R,T))) 



2(tr(W AbPr (R, T))) 
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literature [7,8,10,11]. Our conclusion is therefore that the vacuum type of pure gauge QCD 
with two colors is close to the type I - type II boundary, even if marginally of type II. 

Our result can be further clarified by looking for observable consequences of the QCD 
vacuum being a type I or type II superconductor: that will be the subject of next Section. 

IV. ANALYSIS OF FLUX TUBES INTERACTIONS 

Another direct (but qualitative) method to identify the typology of the vacuum is to 
look at the behavior of two close parallel flux tubes: in a type II superconductor nearby flux 
tubes repel each other, while attraction is expected for a type I superconductor. 

We have looked at the flux tubes of two qq pairs placed with their axes parallel to each 
other (and along the z direction). We have performed two simulations at /3 = 2.6, placing 
each quark at a distance of 16 lattice spacings from the respective antiquark and considering 
two different distances D between the two parallel qq axes, D = 4 and D = 5; the lattice 
chosen is again a 24 x 24 x 32 x 24 and following Section III, we have determined the abelian 
electric field in presence of two parallel Wilson loops Wi(R, T) and W 2 (R,T) at distance 
D = 4aorD = 5a, with R = 16 and T = 6. The quantity we look at is 

(tr (w AbPr (R, T)Wf Pr (R, T) J] $ bPr ) > (tr (wf bPr (R, T)W AbPr (R, T)) tr n i- Pr ) 
El ~ (tr {wf bPr (R, T)W 2 AbPr (R, T)) ) 2 (tr (w^ Pr (R, T)W 2 AbPr (R, T)) ) 

that is the generalization of Eq. (3.1) satisfying the Maxwell equations on the lattice. Our 
statistics consist of about 6 • 10 4 decorrelated configurations for each simulation. 

We focus on the longitudinal E z component in the xz plane reported in Fig. 7 where 
errors on E z are of the order of 5% for almost all points. As in Section III we have used 
cooling for noise reduction: all data showed have been obtained after 6 cooling steps. 

No evident repulsive or attractive behaviour can be appreciated from Fig. 7, however we 
have tried a quantitative analysis of flux tube deflection by measuring the average distance 
between the two tubes and comparing it with the distance D between the two qq axes. The 
average distance has been taken over the central part of the flux tubes, including 9 lattice 
sites for each tube. We have defined the position of the flux tube in two different ways: first 
by the position of the local maximum for E z , secondly by the weighted average position over 
the three lattice sites closest to the qq axis, using E z as a weight. We call the two definitions 
d,M and dw respectively. 

In Table III we report the data obtained for the deflections (d^ — D) and (dw — D) as 
a function of D: a positive/negative value corresponds to a repulsive/attractive behaviour. 
While deflections are nearly compatible with zero at D = 5, some signal appears when the 
flux tubes are closer to each other, at D — 4. Although we have no clear sign of flux tube 
repulsion, we consider this as an important hint in that direction, even more if we consider 
that the superposition of the two flux tubes in the central region should bias our result in 
the opposite direction, leading to negative values of (djw — D) and (dw — D). 

In Section III we concluded that the vacuum is close to the type I - type II boundary, even 
if marginally on the type II side: that would imply a weak repulsive interaction between 
parallel flux tubes. That is consistent with the result of the present Section, i.e. that 
there are signs of weak repulsive interaction as the distance between the two flux tubes is 



9 



decreased. The still large uncertainties as well as the finite lattice spacing place a limit on 
the observable flux tube deflection. While the aim of the study presented in this Section was 
only to look for a possible evident signal of flux tube deformation, more refined investigations 
can be done, including a detailed analysis of Wilson loop interactions and a quantitative 
comparison with the determinations of £ and A, after also taking properly into account the 
quantum fluctuations of the flux tube. That is beyond the purpose of the present study and 
will be the subject of future investigations. 

V. CONCLUSIONS 

The aim of our study was that of understanding which type of dual superconductor 
is realized in the vacuum of SU (2) pure gauge theory. We have followed two different 
strategies: a numerical determination of the parameters £ and A and an analysis of the 
interaction between parallel flux tubes. 

To determine £ we have studied the temporal correlator of an operator fi which creates 
a magnetic monopole and whose vacuum expectation value has been shown to be an order 
parameter for dual superconductivity both in quenched and in full QCD. The greatest diffi- 
culty in our measurement derives from the necessity of isolating an exponentially decaying 
signal from a large background: that is why a very high precision is needed and numerical 
strategies such as analytic link integration have been used. We have determined £ for four 
different values of the inverse coupling (3 and, in some cases, for different lattice volumes 
and for different charges Q of the magnetic monopole. We have explicitely checked that 
our results are independent of the abelian projection used to define /i, that they have no 
significant dependence on the monopole charge Q and that they are not affected by finite 
volume effects within our statistical uncertainties. Data are compatible within errors with 
the correct scaling to the continuum limit, even if a reliable extrapolation to that limit 
cannot still be performed. Results are summarized in Table II and our estimate for the 
correlation length, based on the determination at the largest value of /3, is £ = 0.11 ± 0.02 
fm. We have also repeated our measurement using an alternative way to study the temporal 
correlator [21], obtaining a good agreement. We would like to stress the consistency of our 
results with those reported in Ref . [1 1] , which were obtained through a completely different 
method, consisting in the study of the correlations of monopole currents around the flux 
tube. 

To determine A we employed the usual analysis of the longitudinal component of the 
abelian chromoelectric field inside the flux tube. We performed our measurement at two 
different couplings, (3 = 2.5115 (A = 0.163(7) fm) and (3 = 2.6 (A = 0.160(7) fm), obtaining 
a good agreement with previous literature. 

Our determinations show that £ is smaller than A, even if comparable to it in magnitude: 
this indicates that the vacuum of pure gauge QCD with two colors is close to the type I - 
type II boundary and marginally of type II. This is consistent with our direct investigation 
presented in Section IV, showing some weak signals of repulsive interactions as two parallel 
flux tubes are brought closer to each other. 

One way to improve our results would be to make a more precise determination of £: 
the great noise in the signal obtained for the temporal correlator has been a limitation 
and no significant improvement has been achieved when adopting the alternative approach 
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proposed in Ref. [21]. One reason for the problems encountered can be traced back to the 
small value of £ itself (£ ~ 0.1 fm), which makes the signal to fade away very rapidly after a 
few lattice spacings. One possible strategy to overcome this limitation could be to make use 
of anisotropic lattices with a very small lattice spacing in the temporal direction; we will also 
consider direct determinations of the correlator (p,(t,x)/j,(0,x)), in place of its derivatives 
pit) or p{t), using a technique recently developed for the U(l) gauge theory [27]. Finally, 
a more refined investigation of flux tube interactions could be performed, using also finer 
lattice spacings in order to be more sensitive to small flux tube deflections. 

As an extension of our study, we will repeat in the future the determination of £ also 
around or slightly above the deconfinement critical temperature T c , where £ could be directly 
related to the mass of physical magnetic monopoles: that could be quite relevant from a 
phenomenological point of view, if the hypothesis [28] of light monopole degrees of freedom 
populating the quark gluon plasma slightly above the transition is correct. Of course it 
would be also of fundamental importance to extend our investigation to the case of pure 
gauge theory with 3 colors and eventually to full QCD. 
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TABLES 





to = 1 


t = 2 


to = 3 


t = 4 


to = 5 


12 3 x 16, = 2.4, Q = 2 


0.857 ±0.021 


1.07 ±0.07 


1.02 ±0.20 


4.5 ±3.8 




12 3 x 20, = 2.5115, Q = 2 




1.28 ±0.08 


1.36 ±0.21 


1.8 ± 1.5 




16 3 x 20, = 2.5115, Q = 2 


0.99 ± 0.03 


1.14 ±0.10 


1.28 ±0.25 


1.0 ±0.4 




20 3 x 20, = 2.6, Q = 2 


1.24 ±0.03 


1.33 ±0.05 


1.52 ±0.17 


3.6 ± 1.5 




20 3 x 20, = 2.6, Q = 8 


1.06 ±0.03 


1.32 ±0.08 


1.62 ±0.23 


1.8 ±0.6 


1.2 ± 1.0 


20 3 x 20, = 2.7, Q = 8 


1.34 ±0.02 


1.51 ±0.04 


1.78 ±0.11 


2.4 ±0.4 


1.7 ±0.8 



TABLE I. Dependence of £ on the fit starting point to for our sets of numerical simulations (fit 
according to Eq. 2.10). 






a(0) 




2.4 


0.118 fm 


0.126 ±0.008 fm 


2.5115 


0.083 fm 


0.106 ±0.021 fm 


2.6 


0.062 fm 


0.094 ± 0.011 fm (Q = 2), 0.100 ± 0.014 fm (Q = 8) 


2.7 


0.046 fm 


0.110 ±0.018 fm 



TABLE II. Lattice spacing a and correlation length £ for different values of 0. 



D 


(1m - D 


dw — D 


4 


0.44(17) 


0.017(9) 


5 


0.22(14) 


-0.007(15) 



TABLE III. Average deflection (in lattice spacing units) at = 2.6 for two parallel flux tubes 
as a function of the distance D between the qq axes. See text for the definition of dM and dw- 
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FIGURES 
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FIG. 1. Modified average plaquette as a function of t for a 12 3 x 16 lattice at j3 = 2.4 (top left), 
a 12 3 x 20 lattice at (3 = 2.5115 (top right), a 16 3 x 20 lattice at (3 = 2.5115 (middle left), a 20 3 x 20 
lattice at P = 2.6 with Q = 2 (middle right), a 20 3 x 20 lattice at [5 = 2.6 with Q = 8 (bottom left) 
and a 20 3 x 20 lattice at (3 = 2.7 with Q = 8 (bottom right). In the last case ((3 = 2.7) we have 
actually plotted the quantity log(t 1 / 2 (n — A 1 )) (see Eq. (2.10)) in order to highlight the physical 
signal in the correlator. 
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0.62945 



0.62940 



0.62935 




FIG. 2. Mean modified plaquette in the random gauge and in the Polyakov gauge: an offset as 
been added in order to compare the two sets of data, which lead to compatible correlation lengths. 




FIG. 3. Parameter p measured at j3 = 2.5115 on a 12 3 x 16 lattice. 
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FIG. 4. Profile of E z after 6 cooling steps for a 24 x 24 x 32 x 24 lattice at f3 = 2.6 (left) and 
at (3 = 2.5115 (right). The black line refers to a fit with the function AKq(j). 
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FIG. 5. A as a function of the fit starting point d (left) and as a function of the number of 
cooling steps (right) at = 2.6. 
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FIG. 6. £ and A in fermi units for different values of (3. The two different values shown at 
P = 2.6, which have been slightly split apart for the sake of clarity, correspond to two different 
monopole charges, Q = 8 and Q = 2 respectively. 




FIG. 7. Profile of the two interacting flux tubes on the xz plane placed at different relative 
distances (D = 4 and D = 5): the two qq axes are placed respectively at x = 9 and x = 13 (left), 
and at x = 8 and x = 13 (right). The electric field has been measured after 6 cooling steps. 
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